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I.  INTRODUCTION 

The  Air  Force  Flight  Dynamics  Laboratory  (AFWAL/FI)  has  been  con¬ 
ducting  research  on  the  effects  of  control  system  dynamics  on  longitudinal 
flying  qualities  during  fighter  approach  and  landing.  A  sharp  degradation 
In  flying  qualities  takes  place  during  this  critical  phase  "of  the  landing 
task.  Severe  pilot-induced  oscillations  during  the  flare  have  been 
reported.  The  objective  in  the  program  has  been  to  investigate  pilot- 
induced  oscillations  (PIO)  of  the  NT-33  aircraft  due  to  significant 
control  system  lags,  to  effects  of  wind  shear  and  to  pilot  delays. 

Advance  digital  control  schemes  add  much  greater  flexibility  and  logic 
capabilities  when  compared  to  analog  systems.  However,  this  increase  in 
complexity  of  future  aircraft  flight  control  systems  may  also  add  larger 
control  system  lags.  It  has  been  observed  that  large  control  system 
lags,  high  pilot  gains,  pilot-lag  due  to  neuromuscular  dynamics  and 
aerodynamic  transport  lag  are  all  possible  causes  of  pilot-induced 
oscillation  problems.  These  phenomena  all  require  careful  theoretical 
analysis. 

It  should  be  stressed  that  the  use  of  digital  control  system  is 
now  a  reality  and  Its  effects  on  flying  qualities  of  these  fighter 
aircrafts  need  careful  analysis.  The  variable  stability  NT-33  is 
capable  of  producing  a  wide  range  of  aircraft  and  control  system 
characteristics.  The  main  reason  for  selecting  the  NT-33  aircraft  was 
to  test  the  flying  qualities  of  simulated  YF-12  and  YF-17  aircraft.  The 
simulation  of  the  YF-17  with  the  NT-33  aircraft  has  encountered  serious 
PIO  difficulties  in  flare  whereas  no  such  problems  have  been  reported  for 
YF-17  II] .  Some  detailed  studies  of  PIO  during  the  NT-33  aircraft 
simulation  can  be  found  in  earlier  works  of  USAF/Calspan  {2] .  Calspan 
diagnosed  the  PlO-problem  as  excessive  control  lags.  They  modified  the 
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simulated  control  system  dynamics  to  reduce  the  lag  contribution  to 
longitudinal  dynamics  and  found  that  It  reduced  the  problem.  The  effects 
of  significant  control  dynamics  on  fighter  approach  and  landing  longitudinal 
flying  qualities  were  also  Investigated  In  flight  using  the  USAF/Calspan 
NT-33  aircraft  [3].  Pilot-Induced  oscillations  occurred  during  the 
landing  task.  The  flight  tests  reported  In  [3]  provide  a  data  base  for 
the  development  of  suitable  flying  qualities  requirements  which  are 
applicable  to  aircraft  with  significant  control  system  dynamics. 

The  properties  of  solutions  of  linear  differential  equations  of 
the  retarded  type  with  constant  coefficients  and  constant  time-delays 
for  the  pilot  model  has  been  considered  by  several  authors.  However,  the 
formulation  has  not  been  considered  when  the  coefficients  and  retardations 
in  the  closed-loop  pilot  model  are  variable.  Such  formulation  may  now 
be  justified  when  the  effect  of  wind  shear  and  the  neuromuscular  system 
dynamics  are  Included.  This  extension  of  the  analysis  Is  suggested  by 
the  ••“cent  measurements  that  have  been  cited.  A  generalized  closed- 
loop  nonlinear  pilot  model  for  NT-33  aircraft,  with  variable  retardations 
and  coefficients  Is  considered  herein.  The  theoretical  analysis  Is 
developed  in  the  time  domain  to  analyze  the  pilot-induced  oscillations 
problem  In  the  most  general  format. 

II.  OBJECTIVES 

The  structure  of  the  research  Is  as  follows.  First,  a  formulation 

O 

of  the  closed-loop  NT-33  pilot  model  Is  Introduced.  The  NT-33  air-frame 
dynamic  equations,  linearized  about  the  trim  conditions  and  representing 
the  manual  flare  and  landing  of  the  aircraft,  have  been  used.  The  pilot 
dynamics  are  assumed  to  have  variable  gain  and  variable  retardation, 
possibly  due  to  wind  shear  and  the  neuromuscular  effects.  It  is 
assumed  that  the  pilot  forms  the  closed-loop,  thus  changing  the  overall 
characteristics  of  the  system. 


1 

After  having  introduced  the  required  formulation  of  the  closed- 
loop  NT-33A  pilot  model,  certain  theorems  dealing  with  the  spectral 
factorization  of  entire  functions  of  the  exponential  order  were 
used  to  generate  Lyapunov  functionals.  The  reference  source  for  this 
material  is  the  English  translation  of  the  book  by  Levine  [4J ,  which 
:  gives  a  comprehensive  treatment  of  the  properties  of  the  zeros  of  the 
entire  functions  and  related  topics.  The  spectral  factorization 
theorems  play  a  central  role  in  constructing  the  Lyapunov  functionals. 
Spectral  factorization  is  emphasized  for  the  role  that  these  equations 
play  in  generating  Lyapunov  functionals  for  a  class  of  system  that 
represents  PlO-systems,  rather  than  on  the  mathematical  proofs. 

Convolution  equations  involving  distributions  which  satisfy  assumptions 
made  by  Hale  and  Meyer  J5J  for  the  functional  equations  of  the  delay 
type  are  used  to  describe  the  dynamical  systems  in  analyzing  the  PIO 
system. 

III.  THE  PILOT  MODEL 

Figure  1  represents  a  nonlinear  pilot  model.  The  NT-33  airframe 
dynamic  equations  linearized  about  the  trim  condition  representing  the 
manual  flare  and  landing  of  the  aircraft  have  been  used.  These  equations 
are  the  same  as  those  given  in  USAF/CALSPAN  {2]  and  Smith  [3],  except 
that  only  the  longitudinal  transfer  functions  have  been  derived.  The 
longitudinal  equations  representing  the  open-loop  aircraft  dynamics  about 
the  trim  conditions  during  the  flare  and  landing  maneuvers  are  represented 
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where  symbols  represent  as  follows:  u,  w,  q,  8  are  perturbation  velocity  from 
trim  along  x-body  axis,  perturbation  velocity  from  trim  along  z-body  axis,  body 
axis  pitch  rate  and  pitch  attitude  respectively.  6£S  is  pitch  stick  deflection 
at  grip.  Also,  notationally. 
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Xu.  Zu,  2^  are  body  axis  dimensional  x-force  derivative  and  z-force  deriva¬ 
tives  respectively,  Mu,  Ix,  1^  are  afrcraft  mass,  moment  of  Inertia  about  body 
x-axIs  and  body  y-axls  respectively,  uQ,  wQ,  qQ,  and  0Q  are  trim  values. 

These  equations  Imply  that  the  reference  axis  are  body  axis  and  the  wings 
are  always  level.  For  small  angles,  u0  ®  Vy ,  the  trim  true  airspeed,  and 
°o  *  Wo/VT*  The  variables  u,  w  (a),  8  and  6gS  are  all  Incremental  values  from 

the  reference  trim  conditions. 
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The  transfer  function  for  longitudinal  stick  Input  (°ES)  and  pitch  attitude 
control  output  (6)  can  be  derived  from  the  equations  of  motion  (1)  as 

%s  <S  +  X)<s  *  x1 

A  |  C  /  \  /n\ 


K  <'2  *  2Sp  V  *  "*>>  <’*  *  "ph  *  *  “phi 

where  Tq-j^  represents  airframe  lead  time  constants  csp»  wsp»  and  ?ph,  wph  are 

short  period  damping  ratio,  short  period  undamped  natural  frequency  and  phugoid 
damping  ratio,  phugoid  undamped  natural  frequency  respectively.  p(s)  anda(s) 
are  polynomials  in  s.  The  polynomials  p(s)  ando(s)  are  such  that  the  degree  of 
o(s)  is  assumed  to  exceed  to  that  of  p(s).  An  inspection  of  the  data  base  of 
ref.  [3,  10-12]  suggests  that  a  reasonable  model  for  pilot  dynamics  in  pitch 
tracking  would  be  of  the  form 

L(s)  =  E  a.  (t)  s1  ~sr(t)  +  E  b.  (t)  s*  ;  n  >  m  (3) 

i=o  1  e  1=o  1 

where  a{  (t),  and  b^  (t)  are  bounded  coefficients  and  x(t)  is  a  bounded  time- 
delay,  the  time-delay  T(t)  is  unknown  and  can  be  assumed  due  to  neuromuscular 
effect  of  the  pilot.  A  small  transmission  lag  may  also  be  present.  The  pilot 
dynamics  is  assumed  to  have  variable  gains  a^  (t),  b^  (t)  possibly  due  to  wind 
shear  and  the  neuromuscular  effects.  It  is  assumed  that  the  pilot  forms  the 
closed-loop,  thus  changing  the  overall  characteristics  of  the  system. 

A  closed-loop  analysis  can  be  performed  by  considering  the  pilot  to  be 
controlling  to  some  desired  attitude  which  minimizes  the  pitch  attitude  error  e. 
The  non-linearities  of  the  artificial  feel  system  are  included  in  the  model,  as 
well  as  the  non linearities  in  the  stability  augmentation  systems. 

This  completes  our  formulation  of  the  closed-loop  nonlinear  model  with 
the  pilot  In  the  loop.  The  dynamics  of  the  NT-33  airframe,  the  pilot  and  the 
nonlinearities  of  the  artificial  feel  system  as  well  as  the  nonl Inearl ties  of 
the  stability  augmentation  system  have  all  been  defined.  In  the  next  section, 
we  give  some  notations,  theoretical  backgrounds  and  our  method  of  analysis. 
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IV  APPLICATION  OF  INVARIANCE  PRINCIPLE 

The  'Invariance  principle'  of  J.  P.  LaSalle  played  an  important  role  in 
the  theory  of  abstract  dynamical  systems.  Very  few  practical  applications 
[7,13]  of  'invariance  principle'  have  been  attempted  because  of  its  complexity. 
Nevertheless,  it  is  also  particularly  useful  in  solving  practical  problems. 

By  applying  the  invariance  principle,  we  are  able  to  derive  stability  criteria 
for  the  dynamical  systems  which  are  optimal  in  a  certain  sense.  Thus,  the 
stability  results  can  be  improved  considerable  when  the  study  of  the  system  is 
based  on  LaSalle's  invariance  principle  [6]. 

In  this  section,  we  first  introduce  our  notations.  This  is  followed  by 
some  lemmas  on  spectral  factorization  and  its  application  to  the  invariance 
principle.  Levin  [4]  has  already  given  comprehensive  treatment  of  the  proper¬ 
ties  of  zeros  of  the  entire  functions.  The  spectral  factorization  of  the  entire 
function  plays  an  important  role  in  studying  the  properties  of  solutions  as  we 
shall  see  in  our  analysis  of  the  pilot-induced  oscillation  problem. 

The  following  notation  from  Hale  [8]  will  be  used  in  this  paper:  En  Is 
complex  Euclidean  n-space,  and  for  x  c  En,  |x|  denotes  any  vector  norm.  For 
a  given  r(t)  >_  7  >  0,  C  denotes  the  linear  space  of  continuous  functions  mapping 
the  internal  [-7,0]  into  En  and  for  <p  e  C,  ||$||  =  sup  |  <f>(0)|,  *  7_<e_<0. 
Of  course  C  with  this  norm  is  a  Banach  space.  For  H  >  0,  CH  denotes  the  set  of 
4  in  C  for  which  j|  4>  ||  <  H,  for  any  continuous  function  x(s)  whose  domain  is 
-  7 _<  s  _<ct,  o  >  0,  and  whose  range  is  in  En,  and  for  any  t,  0  _<t  _<a  the 
symbol  xt  will  denote  x^  (0)  ■  x  (t  +  0),  -  7_<e  _<0;  that  is  xt  is  in  CH, 
and  is  that  segment  of  the  function  x(s)  defined  by  letting  s  range  in  the 
interval  t  -  7  _<  s  _<  t. 

Let  G(t,+)  be  a  function  defined  on  (0,«°)  x  Into  En  and  let  x  (t) 
denote  the  right  hand  derivative  of  x(s)  at  s  *  t.  The  system 
x(t)  ■  G(t,  xt),  t  >.  0 

Is  called  a  functional  -  differential  equation  (FOE). 


8 


Definition  -  Let  t0  >,  0,  and  let  $  be  any  given  function  In  C^.  A  function 
x(t0,  $)  (t)  is  said  to  be  solution  of  FDE  with  Initial  function  $  at  t  *  tQ 
If  there  exists  a  number  A  ^  0  such  that 

(a)  for  each  t,  t0_<t  _<t0  +  A,  xt  (tQf  $)  Is  defined  in  C^. 

(b)  xt  (t  ,  ♦)  ■  <t> 

o 

(c)  x(t0,  4»)  (t)  satisfied  the  functional  differential  equation 
(FDE)  for  tQ  _<  t  _<  tQ  +  A. 

To  analyze  the  PlO-problem  the  original  model  in  fig.  1  is  redrawn  as 
e  s  t*e^(t)  P*e2(t)  *  G 


Fig.  2 

Note  that  In  the  block  diagram 
t  *  e1  (t) 

-  (  I  a^t)  bj(t)  6  1  )  *  e^t)  (4) 

Notatlonally,  we  have 

4*  *  e(t)  -  em  (t) 
and 

<5(t)  *  e(t)  -  em(t-T(t)) 

where  *  denotes  the  convolution  operator  and  i*  denotes  the  m-th  derivative 
of  Dirac  delta  function.  Throughout  we  are  assuming  that  the  distribution 
functions  are  measureable  and  have  compact  support  on  [0,t). 
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Next,  we  define  the  convolution  equations, 
e^t)  -  a  *  e2(t) 

e  (t)  ■  £  *  e^(t)  =  l*o*  e2(t)  (5) 

The  Laplace  transform  of  the  distribution  e-j(t)  can  be  represented  by  a 
function 

-st  “>  -st 

<  e^ ,  e  >  =  /  e.|(t)  e  dt 

It  is  assumed,  however,  that  the  distributions  are  of  the  finite  order.  In 
other  words  the  distributions  throughout  have  compact  support.  For  details  of 
such  distribution  functions,  we  refer  to  Schwartz  [9].  Our  objective  here  is 
to  provide  some  background  material  rather  than  vigorous  derivations.  Next, 
we  state  a  Theorem  which  establishes  the  properties  of  solutions  e(t).  This 
result  will  be  used  later  in  analysis  of  the  closed-loop  pilot  model. 


Theorem  1.  The  solutions  of  the  equation 
e(t)  =  l  *  o  *  e2  (t)  =  0 

are  exponentially  stable  provided  the  transcendental  polynomial  t(s)a(s) 
satisfy  the  following  conditions: 

(i)  Re  [l  (s)  o  (s)  |Tjtj  =  o]  <  0 

(11 )  l  (w)  o  (w)  \  0 

Proof  of  Theorem  1.  It  is  sufficient  to  show  that  the  real  parts  of  the 
roots  of  the  transcendental  polynomial 
P(s;x(t))  *  A(s)  •  o(s) 

*  (  Z  a.  (t)  s1  e"T^s  +  Z  b.  (t)s1) 
i*o  1  i*o  1 


•I*2  *  2  ‘sp  "sp  s  ♦  wsp>  ($2  *  2  ‘ph  “ph  s  *  '  0 

are  all  negative  for  t  (t)  >.  0.  It  is  obvious  that  Rea(s)  <0  provided 


(6) 
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0  _5  Csp  _<  wSp  and  0  _<  5ph  _<  wph*  To  show  that  the  real  part  of  the 
transcendental  polynomial  *(s;x(t))  to  be  negative,  we  expand  In  the  form 

£(s;x(t))  =  sn  +  [an  l  (t)  e“T^s  +  bn  l (t)  ]  s  n_1 

+  tan_2  (t)  e"T(t)s  +  bn_2  (t)  ]  s  n'2  +  ... 

+  [a0  (t)  e-T(t)s  +  b0  (t)  ]  (7) 

we  have  assumed  that  m  *  n-1.  If  m  Is  order  less  than  n-1,  we  can  set  the 

coefficients  a„  ,  etc.  zero, 
n- 1 

t(s;x(t))  =  sn  +  Pn_i ( t)  sn_1  +  Pn.2(t)  sn'2  +  ...  +  p0(t)  (8) 

where 

Vi  (t)  =  Vi  (t)  e'T(t>S  +  bn-l  <*> 

P„-2  <*>  ‘  V2  <»)  P'T“)S  ♦  b„-2  <‘> 

•  •  • 

•  •  • 

(t)  -  (t)  e"'(t)s  .  b0  '(t) 

Since,  |  exp  [-t  (t)s ]  |  _<  1  for  all  x(t)  >_  0,  hence,  when  x(t)  >.0  and  s  >  0, 
the  coefficients  p..(t),  i=l,...,  .n  are  bounded.  Let  fT  denote  a  constant  such 
that  p^  =  max  |  p^(t)  |  and  let  D  «  max  [1,  (n+l)p"]  >0.  We  will  now  show  that 

l._<1  _<  n  under  the  assunptions  of  Theorem  1,  all  roots  of  *(s;x(t))  lie  in  the 

left  half  plane.  To  prove  this,  we  consider  two  cases:  (1)  when  |s  |  >_  0  and 
(11)  when  )s|  <  D.  Now  suppose  |s|  >.  D,  then 

I  t(s;x(t))  |  *  |  sn  +  pnl  (t)  sn_1  +....+  p0  (t)  | 

-lsl" 

n  np 

ilsl  t'-TirrTrpl>0  <9> 

The  last  Inequality  follows  from  the  fact  that  |s(  >.  0  ■  (n  ♦  1)  p- >.  1.  Thus, 

In  the  domain  |s|  >  D  and  Re  (s)  >  0  the  pdlynomlal  t(s;x(t))  possesses  no 


root  for  any  bounded  r(t)  >0.  Now  suppose  that  |s|  <  D.  From  condition 
(1)  of  Theorem  1  roots  of  £(s)  o(s)  are  all  In  the  semi-plane  Re  (s)  <0  for 
i(t)  *  0.  Now  when  T(t)  +  0,  the  only  possibility  for  the  characteristic  roots 
to  fall  within  Re  (s)  >  0  Is  that  for  x(t)^0  the  variable  s  runs  along  the 
Imaginary  axis  on  the  s-plane  from  -D  to  D.  But  the  condition  (11)  does  not 
allow  the  roots  to  run  along  the  Imaginary  axis  of  s-plane  and  therefore  under 
our  assumptions  the  characteristic  roots  must  remain  within  the  semi-plane 
Re  (s)  <  0.  This  completes  the  proof. 

For  the  nonlinear  model  In  fig.  2,  the  description  of  the  system  Is 
obtained  as 

l  *  a  *  e2  (t)  +  f  (t,  p*e2  ( t) )  =  u(t);  t  >  0  (10) 

To  analyze  stability  of  the  functional  equation  (10),  we  shall  construct 
a  Lya punov  functional.  The  spectral  factorization  of  the  entire  function  plays 
an  important  role  in  the  construction  of  this  Lyapunov  functional  as  will  be 
seen  later.  We  depend  heavily  upon  the  works  of  Levin  [4]  who  has  shown  that 
spectral  factorization  may  be  applied  to  entire  functions  of  exponential  type. 
We  state  these  results  for  our  convenience. 

Lemma  1.  In  order  that  an  entire  function  F(s)  of  exponential  type  may  be  of 
class  A  It  is  necessary  and  sufficient  that  for  some  fixed  A  >  0  and  for  every 
R  >  A  the  following  Inequality  be  valid: 

R  In  F  (w)  F  (-w)| 

J  - i -  *  ‘ 

A  W£ 

where  ^  Is  a  constant. 

Lemma  2.  For  an  entire  function  F(s)  of  exponential  type  to  have  the 
representation 

F(w)  -  9  (w)  9  (-*) 

where  9  (w)  Is  an  entire  functfion  of  typo  j  (  I  T(t)  |  _<  T)  with  zeros  In 
the  half-plane  Re(s)  >_  0,  If,  and  only  If  F(s)  belongs  to  class  A  and  F(w)  >0. 
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We  will  now  Introduce  our  main  results  to  establish  asymptotic  stability 
of  the  nonlinear  system  (10).  Stability  of  the  largest  invariant  set  is  then 
derived  from  the  application  of  the  'invariance  principle'. 

5.  MAIN  RESULT 

Consider  a  Lyapunov  functional 
V  (•)  -  4Z  q  e2 

t(e9) 

=  /  c  [2  U  *o  *  e?  (t))  (p  *  e2  (t) ) 
t(o)  2  Z 

-U*e2(t))2  ]  dt  (ID 

where  we  have  assumed  that  t(s)  o(s)  and  p(s)  have  no  common  zero.  Let 
the  assumptions  of  lemmas  1  and  2  hold,  then 

def, 

F(s)  *  *(s)  o(s)  p(-s)  +  l(-s)  o(-s)  p(s)  (12) 

has  spectral  factorization  such  that 

F(s)  *  4>(s)  ♦(-s).  03) 

Let  4>(s)  be  the  transform  of  a  distribution  of  order  <  n  and  support  in 
[0,T].  If  the  assumption  of  Theorem  1  holds,  then  F(*)ls.«w>  0.  We  now 
show  that  V(*)  is  positive.  The  state  e2  (t)  «  e2U0»V^  such  that 
at  time  tj  *  t(0),  (e2(t))^  *  0.  And  similarly  at  t2  *  t(e2),  (e2)^  *  e2‘ 
Therefore  t(0)  *  »,  and  t(e2)  *  0.  Hence  the  functional  (11)  can  be 
written  as 

V(0  *  e2  (t)  q  e2  (t) 


-/  (2  U  *  o  *  e2  (t))  (p  *  e2  (t)) 

0 

(♦  *  e2  (t))2  )  dt 
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For  examining  the  Invariance  properties  of  the  nonlinear  system  (15) 
with  u(t)  *  0, 

I  *  a  *  e2  (t)  +  f  (t,  p  *  e2  ( t) )  *  0  (15) 

we  compare  Its  solution  with  the  linear  system 

A  *  a  *  e£  (t)  -  0.  (16) 

Next,  we  show  that  V(*)  >  0  along  the  solution  of  eqn.  (16). 

Obviously 

V(-)  «  /  (♦  *  e2  (t))2  dt  >  0  (17) 

o 

provided  $  *  e2  (t)  0.  We  of  course  assume  that  4>  *  e2  (t)  is 

defined.  Now,  computing  D+  Vt  (?2)  along  the  trajectories  of  eqn.  (16) 
yields: 

D+  Vt  (e2)  «  lim  sup  £  [Vt+h  (?2)  -  Vt  (?2)  ] 

i  _  o 

«  lim  sup  r-  [  /  ($  *  e?  (t))  dt 

h  -  o*  n  t+h  c 

-  I  (♦  *  e?  (t))2  dt  ] 

t  c 

«  -[♦  *  e"2  (t)]2  .  (18) 

Thus  the  solution  of  eqn.  (16)  is  decaying  and  any  two  solutions 
*2t  (to*  V  (t)’  *2t  (to’  no)  of  <16>  sat1sfy  the  estimate, 

II  f2t  (t0.  *„>  -  «2t  (t0.  »0)  ll0_<  A(t.r)||T#  -n0llora(t-tol  tw 

where  notatlonally  \\*\\0  -  "»X_T  *  <0  I W*>l I*  Now  consider 

functional  (17)  and  compute  along  the  trajectories  of  (15) 
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°\('2>(15)  *  F  *  Ve2» 

'  ’iTo*1*1  *  (1Vh<*2(t))  *  W‘2<‘>» 

♦  [*t+h(ij(t))  -  »t  (e2(t))]  ) 

<  llmsup  1  mvt,h(.2(t))  ■  *tth(?2(t))||  ♦  d\(?2)„6) 

<-[**e2(t)]2  +  C  ||f(t,p  *  e2(t)|| 

<  -  t4>  *  e2(t)]2  +  Cg(t,||e2  (t)||Q) 

<  -  W(t,e2(t0^0)(t))  (20) 

We  have  assumed  that  the  nonlinear  system  admits  unique  solution  and  the 
nonlinearity  is  such  the  ||f(t,p  *  e2(t))| |»<g(t,| |e2<.(t)|(0)  where 

g(t,u)  Is  nondecreasing  In  u  for  tfeR+  and  W  Is  some  non  negative  function. 

It  Is  of  course,  possible  to  obtain  the  function  W  (•)  If  the  distributions 
Involved  are  restricted  further.  However,  a  few  examples  will  be  presented 
to  Illustrate  the  application  of  the  method  and  Indicate  the  type  of  functionals 
that  are  obtained.  To  examine  the  Invariance  properties  of  the  nonlinear 
system,  we  first  note  that  the  solution  e2  (tQ,  rpQ ) ( t )  Is  bounded  for  all 
t  >  tQ  because  of  the  existence  of  a  Lypunov  functional.  If  the  limit  point 

ft  of  e2  (tQ,  4»0)  (t)  does  not  exist  on  the  boundary  G,  G  Is  an  open  set  In  C, 

then  by  lemma  4.8  In  [6],  e2  (t0»  4*0)(t)  Is  pre-compact.  Now  define  the 

largest  Invariant  set  ft  C  E  *  (e2(t0»4»0)(t)  :  W(t,e2(t0^0)(t)  *  0  ) 
then 

for  soam  C  and  all  t  >  o. 
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For  our  assertion,  we  are  making  use  of  the  Theorem  4.7,  P.78  In 
[6].  Thus,  these  results  are  combined  Into  the  following  theorem. 

Theorem  2.  If  there  exists  $(s)  whose  distribution  Is  of  order  and 
compact  support  In  [0,T]  such  that  F(s)  defined  In  eqn.  12  has  the 
spectral  factorization  such  that 

1)  F(s)  ■  *(s)  *(-s); 

11)  Conditions  of  Theorem  1,  and  lemmas  1  and  2  hold; 

111)  ||f(t,  p  *  e2(t))  ||<g(t,  ||e2t(t)||e 

where  g(t,u)  Is  non  decreasing  in  u  for  tcR+ 

Then  there  exists  a  Lyapunov  functional  Vt( • )  on  G,  for  system  (10)  with 

u(t)  *  0.  Furthermore  If 

dtf. 

1  v)  8CE  «  <e2U0.  *o)(t)  :  -  0  } 

where  W  is  some  non-negative  function  then 

*2<V  VW  V‘1(c) 


for  some  c  and  t  >  o. 


I 


*wzpT>y.'  i  .wfM,  wt«V)iw^r»MWW>w-» r*«- - 
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